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Abstract: 



Derived Meteorological Products (DMPs, including potential temperature ( 6 ), potential vorticity, 
equivalent latitude (EqL), horizontal winds and tropopause locations) have been produced for the 
locations and times of measurements by several solar occultation (SO) instruments and the Aura 
Microwave Limb Sounder (MLS). DMPs are calculated from several meteorological analyses for 
the Atmospheric Chemistry Experiment-Fourier Transform Spectrometer, Stratospheric Aerosol 
and Gas Experiment II and III, Halogen Occultation Experiment, and Polar Ozone and Aerosol 
Measurement II and III SO instruments and MLS. Time-series comparisons of MLS version 1.5 
and SO data using DMPs show good qualitative agreement in time evolution of O 3 , N 2 O, H 2 O, 
CO, HNO 3 , HC1 and temperature; quantitative agreement is good in most cases. EqL-coordinate 
comparisons of MLS version 2.2 and SO data show good quantitative agreement throughout the 
stratosphere for most of these species, with significant biases for a few species in localized 
regions. Comparisons in EqL coordinates of MLS and SO data, and of SO data with 
geographically coincident MLS data provide insight into where and how sampling effects are 
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important in interpretation of the sparse SO data, thus assisting in fully utilizing the SO data in 
scientific studies and comparisons with other sparse datasets. The DMPs are valuable for 
scientific studies and to facilitate validation of non-coincident measurements. 


Popular Summary 

This paper is associated with comparisons of sparse space-based observations of the structure of 
Earth’s middle atmosphere. 

Solar occultation techniques provide a method of obtaining highly accurate profile information 
on temperature and trace gases in the middle atmosphere. The method involves looking at the 
sun on layers of near-horizontal paths through the atmosphere, including one view of the sun 
above the atmosphere, which provides a means of calibration for the observing instrument. This 
opportunity to obtain the unimpeded view of the sun on each profile is one of the reasons for the 
high accuracy of the technique, because this knowledge can be used to continually correct for 
changes in solar output at the wavelengths observed and for any drifts in instrument 
performance. 

A disadvantage of solar occultation techniques is that they tend to provide very sparse data, with 
only several tens of profiles per day. This means that cross-validation of data from platforms in 
different orbits is complicated because of lack of spatio-temporal coincidence of these 
observations. Sophisticated interpolation techniques are necessary to obtain meaningful 
comparisons between instruments. They are also needed to complete chemical budgets, if (for 
example) some constituents are observed from one platform while other related gases are 
observed from a different platform. 

One technique for validation is to use meteorological analyses to establish “regimes” where the 
constituent measurements are likely to be related to each other. This approach relies on some 
reasonable assumptions about correlations among different gases, which are borne out by the fact 
that long-lived trace gases, whose distributions are primarily determined by local transport, 
correlate well when the correlation is performed on “material” surfaces. These surfaces, which 
are potential temperature in the vertical and potential vorticity in the horizontal, allow the trace 
gas distributions to be compiled effectively as latitude-height distributions in potential- 
temperature/potential-vorticity space. The potential vorticity component is converted to 
“equivalent latitude” which is that latitude that each potential vorticity contour would enclose if 
it were distributed symmetrically around the globe. 

This paper presents a thorough validation of measurements from a number of solar occultation 
instruments, using “derived meteorological products” (the potential temperature and potential 
vorticity) from Goddard Earth Observing System, Version 5 (GEOS-5) and other meteorological 
analyses. Data from the limb-sounding Microwave Limb Sounder (MLS) on NASA’s EOS-Aura 
platform are also used. The results reveal generally very good agreement among the MLS data 
and SO datasets, with some anomalies for specific gases in certain regions of the atmosphere. 
These results attest not only to the quality of the observations, but also to the physical robustness 
of the GEOS-5 and other meteorological analyses. 
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Abstract. Derived Meteorological Products (DMPs, including potential tem- 
perature (0), potential vorticity, equivalent latitude (EqL), horizontal winds and 
tropopause locations) have been produced for the locations and times of measure- 
ments by several solar occultation (SO) instruments and the Aura Microwave Limb 
Sounder (MLS). DMPs are calculated from several meteorological analyses for 
the Atmospheric Chemistry Experiment-Fourier Transform Spectrometer, Strato- 
spheric Aerosol and Gas Experiment II and HI, Halogen Occultation Experiment, 
and Polar Ozone and Aerosol Measurement II and III SO instruments and MLS. 
Time-series comparisons of MLS version 1.5 and SO data using DMPs show 
good qualitative agreement in time evolution of O 3 , N 2 O, H 2 O, CO, HNO 3 , HC1 
and temperature; quantitative agreement is good in most cases. EqL-coordinate 
comparisons of MLS version 2.2 and SO data show good quantitative agreement 
throughout the stratosphere for most of these species, with significant biases for 
a few species in localized regions. Comparisons in EqL coordinates of MLS 
and SO data, and of SO data with geographically coincident MLS data provide 
insight into where and how sampling effects are important in interpretation of the 
sparse SO data, thus assisting in fully utilizing the SO data in scientific studies 
and comparisons with other sparse datasets. The DMPs are valuable for scientific 
studies and to facilitate validation of non-coincident measurements. 
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1. Introduction 

The characterization of atmospheric observations by air 
mass properties, such as location with respect to the strato- 
spheric polar vortex or the tropopause, is an invaluable tool 
in numerous analyses, both for research studies and valida- 
tion. It is especially valuable, and has been widely used, 
in studies using solar occultation (SO) satellite data, which 
consist of no more than 15 profiles per day at each of two 
latitudes. Bute hart and Remsberg [1986] and halt et al 
[1990] mapped atmospheric trace gases with respect to po- 
tential vorticity (PV) and equivalent latitude (EqL, the lat- 
itude that would enclose the same area between it and the 
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pole as a given PV contour, Butchart and Remsberg [1986]). 
Manney et al [1999] used PV and EqL with satellite data 
to enable detailed study of polar vortex dynamics and as- 
sociated transport in ATMOS (Atmospheric Trace Molecule 
Spectroscopy) data from the ATLAS (Atmospheric Labora- 
tory for Applications and Science) space-shuttle missions; 
a number of other studies of ATMOS, Halogen Occupa- 
tion Experiment (HALOE), Polar Ozone and Aerosol Mea- 
surement (POAM) II and III, and Stratospheric Aerosol 
and Gas Experiment (SAGE) II and III data [Schoeberl 
et al , 1995; Randall et al , 2005, and references therein] 
have used EqL or PV to help realize the full “condition- 
space” coverage of the sparse SO datasets. Other stud- 
ies of limb-sounding datasets, such as the Microwave Limb 
Sounder (MLS) instruments on the Upper Atmosphere Re- 
search Satellite (UARS) and Earth Observing System (EOS) 
Aura missions, have used PV or EqL mapping to study po- 
lar vortex dynamics and trace gas evolution [Manney et al . , 
1995a, 2005a, and references therein]. 

Air mass characterization also facilitates comparisons us- 
ing measurements that are geographically sparse and may 
not fulfill traditional coincidence criteria based on close 
matching of time and location. Santee et al. [2007a, b, this 
issue] use EqL as a coordinate to compare Aura MLS data 
with Upper Atmosphere Research Satellite (UARS) MLS 
measurements taken in the 1990s. EqL and PV mapping 
have been used in studies of aircraft and ground-based ob- 
servations [e.g., Lait et al. , 1990; Redaelli et al ., 1994]. 
Manney et al [2001] used EqL mapping, standard geo- 
graphical coincidence criteria augmented by PV matching, 
and trajectory histories to compare ozone from seven instru- 
ments during the November 1994 period of the ATLAS-3 
mission, including four SO instruments (ATMOS, HALOE, 
SAGE II and POAM II). PV or EqL can be used to augment 
or replace conventional coincidence criteria [e.g., Michelsen 
et al. , 2002; Lumpe et al. , 2003; Thomason et al ., 2006]. 

Several SO datasets are available for validation of 
Aura measurements, including HALOE, SAGE II and III, 
POAM III, and the Atmospheric Chemistry Experiment 
Fourier Transform Spectroscopy (ACE-FTS) instrument. To 
facilitate non-coincident validation and intercomparison of 
measurements sorted by air-mass characteristics, and for use 
in research studies combining Aura with SO datasets, a set 
of “derived meteorological products'’ (DMPs) has been cal- 
culated for these SO datasets; the DMPs consist of fields 
derived from the meteorological analyses’ winds and tem- 
peratures, such as PV and EqL, and fields from the me- 
teorological analyses, all interpolated to the locations and 
time of the satellite observations. DMPs have also been cal- 
culated for Aura MLS version 1.5 (vl.5) and version 2.2 
(v2.2) datasets; as well as contributing to validation and sci- 


ence studies, these are used in producing plots for routine 
inspection of MLS data, and daily EqL/potential tempera- 
ture (0) cross-sections that are posted on the MLS website 
[http://mls.jpl.nasa.gov]. In the following, we document the 
DMPs for the SO instruments and MLS and use them for 
comparisons of SO and MLS data. We explore effects of 
the satellites’ different sampling and extend more traditional 
intercomparisons for validation. DMPs for MLS and/or SO 
instruments are also used in validation and comparisons in 
other papers in this issue [Pawson et al ., 2007; Santee et al ., 
2007a, b]. 

2. Dataset Descriptions 

2.1. Solar Occultation Datasets 

2.1.1. ACE-FTS SCISAT-1, otherwise known as the 
Atmospheric Chemistry Experiment (ACE) [Bernath et al, 
2005] was launched August 2003. The primary instru- 
ment is the ACE-FTS (hereinafter referred to as ACE), 
a Fourier transform spectrometer featuring high resolution 
(0.02 cm -1 , corresponding to a ±25 cm maximum optical 
path difference) and broad spectral coverage in the infrared 
(750-4400 cm -1 ). ACE works primarily in the solar occulta- 
tion mode, collecting atmospheric limb measurements using 
the sun as a radiation source. Version 2.2 of the ACE re- 
trievals [Boone et al , 2005] is used here, except for O 3 , for 
which the ACE product known as “version 2.2 ozone up- 
date” is used. Early validation efforts with the ACE data 
identified a roughly 10% low bias for altitudes near the O 3 
concentration peak when compared to other satellite mea- 
surements [Walker et al , 2005; Fussen et al , 2005; Petelina 
et al. , 2005]. Version 2.2 processing uses microwindows 
in two spectroscopic regions: 1000-1150 cm -1 and 1800- 
2150 cm“ l . Version 2.2 O 3 update uses microwindows in 
the 950-1150 cm -1 range and, in preliminary comparisons, 
exhibits improved agreement with other datasets near the O 3 
concentration peak. ACE vertical resolution is ^3-4 km. 
Latitudes of measurements vary over an annual cycle with 
coverage as high as ±85° and an emphasis on the polar 
regions in winter and spring. Separate files with ACE ge- 
olocation information, primarily latitude and longitude as a 
function of altitude, are provided for each occultation; for 
occultations with missing geolocation files, the geolocation 
information is taken from the headers of the data files, which 
give 30-km tangent point latitude and longitude values; most 
differences are small, but can be up to ~3° latitude and ~ 10 ° 
longitude for brief periods, depending on the viewing geom- 
etry. 

2.1.2. HALOE HALOE [Russell et al , 1993] was 
operational on the Upper Atmosphere Research Satel- 
lite (UARS) from October 1991 through November 2005. 
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HALOE observations take approximately one month to 
cover the full range of latitudes sampled (ranging from ±80° 
to ±50°, depending on season, reaching highest (lowest) lat- 
itudes in summer (winter)). Estimates of uncertainties for 
the profiles of the retrieved HALOE parameters from its first 
public release dataset (Version 17) are provided in Journal 
of Geophysical Research , Atmospheres , 101, D6 - a special 
UARS Validation Issue published in 1996. The data used 
here are Version 19. Some updated uncertainty estimates 
for Version 19 are available for O 3 [Randall et al. , 2003], 
H 2 O [Kley et al , 2000], and temperature Remsberg et al 
[2002]. Vertical resolution is ~ 2 km for O 3 and H 2 O, ~3- 

4 km for temperature, and ~4 km for HC1. Temperatures 
in the HALOE files below 35 km are from the NCEP/CPC 
meteorological analyses. HALOE latitude and longitude as 
a function of height are provided on the same 0.3 -km grid 
used for the temperature files. 

2.1.3. POAM POAM II [Glaccum et al. , 1996] and 
POAM III [Lucke etal. , 1999] were visible/near-infrared so- 
lar occultation instruments that typically made 14-15 mea- 
surements per day in each hemisphere around a circle of lati- 
tude with a longitudinal spacing of about 25°. The latitudinal 
coverage was identical each year, slowly varying between 
54°N - 71°N and 65°S - 88 °S. POAM II obtained data from 
October 1993 until November 1996, when the host satellite 
failed; POAM III obtained data from late April 1998 through 
early December 2005. POAM II provided measurements of 
O 3 , aerosol extinction, and NO 2 ; POAM III provided H 2 O in 
addition to these. The POAM III Version 4 ozone retrievals 
differ little from the Version 3 retrievals described by Lumpe 
et al. [2002], and validated by Lumpe et al. [2003] and Ran- 
dall et al. [2003]. At 15 km and above, the O 3 retrievals have 
a vertical resolution of ^1 km and an estimated precision of 
5% [Lumpe et al ., 2002]. The H 2 O retrievals extend from 

5 to 50 km with 5-7% precision and a vertical resolution 
ranging from 1 km in the lower stratosphere to 3 km in the 
upper stratosphere. The H 2 O retrievals have been validated 
by Lumpe et al. [2006]. Geolocation information (latitude, 
longitude and line-of-sight (LOS) angle) for POAM III were 
calculated on an 8 -km grid, and interpolated linearly to the 
POAM 1 -km measurement grid. 

2.1.4. SAGE II SAGE II 

[http://science.hq.nasa.gov/missions/satellite_45.htm] 

used radiances at 600 nm to derive O 3 and at 940 nm to de- 
rive H 2 O; it took measurements from October 1984 through 
August 2005. The instrument and earlier versions of the 
retrieval algorithm, as well as O 3 validation, are discussed 
by Chu etal. [1989], Cunnold et al [1989], and McCormick 
et al. [1989]. SAGE II data used here are Version 6.2. O 3 
data have km or less vertical resolution, and H 2 O data 
no better than ~1 km. Validation of v6.2 O 3 and H 2 O is 


discussed by Wang et al. [2002] and Taha et al. [2004], 
respectively. The precision of SAGE II O 3 is estimated to 
be ~2% [Borchi and Pommereau, 2006] and no credible 
estimates of the precision of H 2 O measurements exist. 
The SAGE II coverage follows a pattern similar to that of 
HALOE. Geolocation information, including line-of-sight 
(LOS) angle, was provided on the measurement grid in 
separate files. 

2.1.5. SAGE III The SAGE III instrument 
[http://science.hq.nasa.gov/missions/satelliteiLhtm] used 
radiances from several channels in the 570 to 600 nm region 
to derive O 3 and at several more in the 925 to 960 nm 
region to derive H 2 O; it took measurements from May 2002 
through December 2005. The instrument and the retrieval 
algorithm, as well as O 3 validation, are discussed by 
Mauldin et al. [1985], McCormick et al. [2002], and Wang 
et al. [2006]. SAGE III data used here are Version 3. O 3 and 
H 2 O data have a vertical resolution of ~1 km. The precision 
of SAGE III O 3 has not been objectively estimated but 
should be comparable to SAGE II. The SAGE III coverage 
from the sun synchronous orbit is confined to mid to high 
latitudes with sunrise events in the southern hemisphere 
and sunset events in the north. Geolocation information 
is included in the SAGE III data files every 10 km; this is 
interpolated linearly to the 0.5-km data grid. 

2.2. Aura MLS Dataset 

MLS measures millimeter- and submillimeter- 
wavelength thermal emission from the limb of Earth’s 
atmosphere. Detailed information on the measurement 
technique and the MLS instrument on the EOS Aura 
satellite is given by Waters et al. [2006]. The Aura MLS 
fields-of-view point in the direction of orbital motion and 
vertically scan the limb in the orbit plane, leading to data 
coverage from 82°S to 82°N latitude on every orbit. Vertical 
profiles are measured every 165 km along the suborbital 
track and have a horizontal resolution of ~200-300km 
along-track and ^3-9 km across-track. Vertical resolution 
of the Aura MLS data is ~ 3-4 km in the lower and middle 
stratosphere, depending on the product [Froidevaux et al . , 
2006; Livesey etal., 2005, available from the MLS web site, 
http://mls.jpl.nasa.gov] . 

Numerous species are retrieved from MLS observations; 
examples using DMPs to compare many of these with the SO 
datasets are presented here. Further details on these species 
are given in papers in this issue: Lambert et al. [2007] (N 2 O 
and H 2 O), Froidevaux et al. [2007b], Jiang et al. [2007], 
and references therein (O 3 ), Froidevaux et al. [2007a] (HC1), 
Santee et al. [2007b] (HNO 3 ), Pumphrey et al. [2007] (CO), 
and Schwartz et al. [2007] (temperature). These papers val- 
idate version 2.2 (v2.2) MLS data and discuss changes from 
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version 1.5 (vl.5); initial validation comparisons of vl.5 
MLS data are given by Froidevaux et al. [2006] and Bar- 
rett et al. [2006]. Reprocessing with MLS v2.2 is ongoing 
and will take over a year to complete; a small subset of the 
MLS data spaced throughout the mission have been repro- 
cessed in v2.2. Several of the comparisons shown here focus 
on reprocessed periods, but we also show some comparisons 
with vl.5 for examples involving timeseries for which v2.2 
data are not yet available. 

23. Meteorological Analyses 

DMPs are currently calculated from the Met Office 
(MetO) dataset for each of the instruments, from the NASA 
Global Modeling and Assimilation Office’s (GMAO) God- 
dard Earth Observing System (GEOS) datasets (4 and/or 5) 
for MLS and ACE, and from the NCEP/CPC (National Cen- 
ters for Environmental Prediction/Climate Prediction Cen- 
ter) analyses for SAGE II. A brief description of these 
datasets follows; further information is given by Manney 
et al. [2005b, and references therein], which also pro- 
vide comparisons between these and other meteorological 
datasets. Pawson et al. [2007, this issue] also show some 
comparisons for the Aura mission time period of GEOS-4 
and 5 with other analyses. 

The MetO data through 12 March 2006 are from the 
stratosphere-troposphere (STT) data assimilation system 
first developed for the Upper Atmosphere Research Satel- 
lite (UARS) project [Swinbank and O’Neill, 1994], and have 
been produced since October 1991. The assimilation used 
an analysis-correction scheme as described by Lorenc et al. 
[1991] until late 2000, when a three-dimensional variational 
(3D-Var) scheme was implemented [Lorenc et al., 2000]. In 
late 2003, a new dynamical core [Davies et al., 2005] was 
implemented in the assimilation system [Swinbank et al, 
2002, 2004]. The MetO-STT data (three-dimensional winds, 
temperature, and geopotential height) are supplied once- 
daily at 12UT on a 2.5° latitude by 3.75° longitude grid, and 
at UARS pressure levels (6 levels per decade in pressure) 
between 1000 and 0.3 hPa (0.1 hPa after late 2003). Af- 
ter 12 March 2006, the stratospheric analyses are provided 
from the same numerical weather prediction (NWP) model 
system as operational forecasts from the Met Office [D. Wal- 
ters, et al, “Enhancing Vertical and Horizontal Resolution 
in the Met Office Global NWP (Unified) Model, in prepara- 
tion]; the same fields are provided, but on a 0.375° latitude 
by 0.5625° longitude grid, at 27 levels from 1000 to 0.4 hPa. 
When there is a need to distinguish, the recent MetO analy- 
ses are referred to as MetO-NWP and the original ones de- 
veloped for UARS as MetO-STT. DMPs are calculated for 
all instruments from the MetO data. 

The GEOS-4 analyses are described by Bloom et al. 


[2005]; a Physical Space Statistical Analysis Scheme is 
used. The GEOS-4 data used here are provided on 55 hy- 
brid (a/pressure) model levels from the surface to 0.01 hPa. 
The horizontal grid is 1.0° latitude by 1.25° longitude, six- 
hourly average fields are provided centered at 0, 6, 12 
and 18 UT. Besides the standard meteorological variables, 
GEOS-4 products include an extensive set of fields from the 
model and assimilation system, including PV calculated in- 
ternally in the model. DMPs for ACE and for MLS vl .5 data 
are calculated from GEOS-4. 

GEOS-5 analyses [Reinecker et al., 2007] have been pro- 
duced for the full period of the Aura mission, and will even- 
tually replace GEOS-4. GEOS-5 uses the Gridpoint Statis- 
tical Analysis method of Wu et al. [2002], a 3D-Var system, 
with a six-hour analyses window. Analyses are produced for 
surface pressure, temperature, winds, moisture and ozone. 
Along with operational meteorological products, the AIRS 
data stream from EOS-Aqua was assimilated, as described 
in Stajner et al. [2007]. The interface between the obser- 
vations and the GCM is performed using the incremental 
analysis update (IAU) approach [Bloom et al., 1996], which 
avoids shocking the model, thus producing smoother analy- 
ses. An evaluation of some aspects of GEOS-5 relevant to 
the Aura mission is presented in Pawson et al. [2007, this 
issue]. GEOS-5 analyses are provided on 72 model levels 
from the surface to 0.01 hPa, and a 0.5° latitude by 2/3° 
longitude grid. DMPs for ACE and for MLS v2.2 data are 
calculated from GEOS-5. 

NCEP/CPC analyses (used for SAGE II DMPs) are from 
an objective analysis at levels above 10 hPa (above 100 hPa 
prior to April 2001) [Finger et al., 1965, 1993; Gelman 
et al., 1986, 1994]; these analyses have been available since 
June 1979. Analyses at and below 100 hPa are from the tro- 
pospheric analysis and forecast cycle [e.g., Derber and Wu, 
1998; McNally et al, 2000]. The NCEP data are provided 
once a day at 12UT on a 65x65 point polar stereographic 
grid for each hemisphere; the fields used here are interpo- 
lated to a 2.5° x 5° latitude/longitude grid. Only temperature 
and geopotential height are provided in the stratosphere, so 
horizontal winds are calculated from the NCEP geopoten- 
tial heights using a form of the primitive equations that ne- 
glects the vertical advection and time tendency terms [Ran- 
del, 1987; Newman et al, 1989]. 

3. DMP Field Description 

Table 1 lists the DMPs calculated for the SO instru- 
ments; these DMPs are produced on the vertical grids used 
for the SO instruments’ data. Table 2 lists the DMPs pro- 
vided for the MLS instrument. MLS records over 100 times 
more profiles per day than the SO instruments, so calculat- 
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Table 1 . Derived Meteorological Product (DMP) Fields and Units for Solar Occultation Instruments 


Field 

Units 

Description 

Geolocation (from Instrument Geometry) 

Alt 

km 

(2D 1 ) Altitudes 

Lat 

deg 

(2D) Latitudes as a function of Altitude 

Lon 

deg 

(2D) Longitudes as a function of Altitude 

Sun Dir 2 

deg cw from N 

(2D) Line-of-Sight angle (LOS) 

Interpolated from Meteorological Data 

Temperature 

K 

(2D) Temperature from meteorological data 

Geop Hgt 

m 

(2D) Geopotential Height 

Zonal Wind 

m/s 

(2D) Zonal Wind 

Merid Wind 

m/s 

(2D) Meridional Wind 

Calculated from Meteorological Data 

0 

K 

(2D) Potential Temperature from met data 

PV 3 

10~ 4 Km 2 kg -1 s- 1 

(2D) Potential Vorticity 

Scaled PV 

10 -4 s~ l 

(2D) Scaled PV, in “vorticity units” 4 

EqL 

deg 

(2D) Equivalent Latitude 

Hor PV Grad 

- 

(2D) Normalized horizontal PV gradient 

Hor T Grad 

K/km 

(2D) Horizontal temperature gradient 

LOS T Grad 2 

K/km 

(2D) Temperature gradient along LOS 

LOS PV Grad 2 

(10“ 4 Km 2 kg- 1 s- 1 )/km 

(2D) PV gradient along LOS 

EqL - VEC 

deg 

(2D) Distance (EqL) from vortex edge center 

EqL - VEI 

deg 

(2D) Distance (EqL) from inner vortex edge 

EqL - VEO 

deg 

(2D) Distance (EqL) from outer vortex edge 

Dyn Tropopause 

km 

(ID 1 ) Dynamical tropopause altitude 5 

TG Tropopause 

km 

(ID) WMO tropopause altitude 


1 2D indicates profile information, ID a single value for each occultation 

2 Not Available for ACE or HALOE DMPs 

3 Interpolated directly from provided dataset for DMPs derived from GEOS -4 

4 Dunkerton and Delisi [1986]; Manney etal. [1994b] 

5 3.5 PVU joined to 380 K 0 in tropics, see text 
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ing the DMPs is computationally intensive. To make the 
calculations feasible, and since the MLS positions are not 
altitude-dependent, MLS DMPs are calculated and output 
on standard 0 or pressure levels, allowing the calculations to 
be done in advance on the gridded meteorological analysis 
fields once for each analysis time, and then interpolated to 
the MLS times and positions. File formats and access to the 
DMPs are described in the Appendix; the DMP calculations 
are described in more detail below. 

3.1* Description of Calculations and Interpolations 

Interpolation of fields provided in the meteorologi- 
cal analyses (horizontal winds, temperature, geopotential 
height, PV from GEOS-4 and GEOS-5 analyses) is done lin- 
early in time (between six-hourly fields for the GEOS anal- 
yses, daily fields for the MetO analyses) and bi-linearly in 
latitude and longitude. Vertical interpolations are linear in 
log(0) for PV or log(pressure) for the other products. EqL is 
calculated on isentropic surfaces and interpolated linearly in 
log(0); for the high resolution MetO-NWP and GEOS prod- 
ucts, the computationally intensive EqL calculation is done 
once and EqL saved on standard isentropic surfaces; these 
files are read in and interpolated as for PV. For the MetO and 
NCEP DMPs, PV is calculated as described by Manney et al 
[1996b], based on the algorithm of Newman et al [1989]; 
PV is provided from the assimilation model in the GEOS 
datasets, and this is used in the DMP files, as it is more 
fully consistent with the analyzed fields than an offline cal- 
culation. Scaled PV (sPV) is in “vorticity units” [10 -4 s _1 , 
Bunker ton and Delisi, 1986; Manney et al., 1994b], giving a 
similar range of values at levels throughout the stratosphere. 

Horizontal PV gradients for the SO instruments are cal- 
culated using PV interpolated to the isentropic level of each 
measurement in the profile, and on the standard 0 surfaces 
for MLS. The PV gradients provided in the DMP files are 
normalized to the hemispheric average on isentropic sur- 
faces; the calculation is also strongly dependent on the reso- 
lution of the meteorological analysis used; thus, the magni- 
tudes of gradients calculated from different meteorological 
analyses are not directly comparable. Horizontal tempera- 
ture gradients for the SO instruments are calculated on the 
pressure level of each measurement. LOS PV and tempera- 
ture gradients are provided for several of the SO instruments 
for which the LOS angle information is available. These 
are potentially useful for validation and data quality stud- 
ies in assessing the homogeneity of atmospheric conditions 
along the LOS. As with the horizontal gradients, the LOS 
PV (temperature) gradients are calculated on the 0 (pressure) 
level of each measurement. 

The vortex edge “center” position is defined as the EqL of 
the maximum of the windspeed times the PV gradient. This 


is very similar, but not identical, to the widely used criterion 
described by Nash et al. [1996]. The “inner” and “outer” 
vortex edges are EqL on vortex and extravortex sides, re- 
spectively, where that vortex definition function changes cur- 
vature, as described for PV gradient by Nash et al [1996]. 
The vortex is considered to be undefined if 0 < 345 K, or 
if any of the following apply: The EqL value at the vor- 
tex edge center is greater than 80° (vortex too small, e.g., 
Manney et al. [1994a]), the windspeed is less than 15.2 m/s 
(polar night jet is too weak, e.g., Nash et al. [1996]), or the 
normalized PV gradient is less than 1.1 (PV gradient is not 
significantly above average). Figure 1 shows the windspeed, 
PV gradient, and calculated vortex edge center position from 
MLS GEOS-4 DMPs during northern hemisphere (NH) win- 
ter and spring (SH winter cases, not shown, are less ambigu- 
ous since vortex is strong and simply-defined). Every day is 
included in the monthly averages for PV gradient and wind 
speed, but only those days on which a vortex edge is defined 
for the position relative to the vortex edge; thus, in some 
cases (e.g., SH lower stratosphere in March, when vortex is 
just starting to develop and does so substantially over the 
month) there is a vortex edge defined when the plotted PV 
gradients and windspeeds do not show an obvious transport 
barrier. The inclusion of windspeed in the definition reduces 
the likelihood of spurious maxima in the PV gradient at high 
EqL being identified as the vortex edge, and the use of the 
combined function provides a means for choosing which of 
double peaks in windspeed and PV gradient [common in fall, 
e.g., Manney et al. , 2002] are selected. Note that at the low- 
est levels in January, the top of the upper tropospheric sub- 
tropical jet is mis-identified as the vortex edge; this is a com- 
mon occurrence, but setting the 0 limit for the vortex to be 
defined higher would eliminate much of the SH sub vortex re- 
gion in winter and spring. In the upper stratosphere, and near 
the stratopause, the jet/P V gradient structure is much more 
complex, and thus mis-identification (e.g., January in SH) is 
common; in fact, there is often not a single most appropriate 
definition of the vortex edge here. Part of the ambiguity is in 
the PV field, because of the decreased static stability in the 
region surrounding the stratopause; far above the cold pool, 
this layer is often nearly isothermal. 

Automated vortex edge definition is most robust under 
conditions for which the vortex is well, and simply, defined, 
i.e., when there is a single region of strong PV gradients as- 
sociated with the polar night jet. This is the case in the mid- 
dle to lower stratosphere during SH winter and many NH 
winters; in these cases the definition used here agrees closely 
with that of Nash et al [1996] and with other methods of 
determining the vortex edge (including simply using a rep- 
resentative sPV contour). In the upper stratosphere, vortex- 
edge location is more ambiguous and dependent on the spe- 
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Table 2. Derived Meteorological Product (DMP) Fields for MLS 


Field 

Units 

Description 

On Standard Pressure Surfaces 

Pressure 

hPa 

Standard Pressure Level Values 

e 

K 

Potential Temperature from met data 

Geop Hgt 

m 

Geopotential Height 

On Standard Potential Temperature Surfaces 

0 

K 

Standard 0 level values 

Pressure 

hPa 

Pressure on 0 surfaces 

Zonal Wind 

m/s 

Zonal Wind 

Merid Wind 

m/s 

Meridional Wind 

PV 

1q- 4 ]^ m 2 ^g-i s -i 

Potential Vorticity 

EqL 

deg 

Equivalent Latitude 

Hor PV Grad 

- 

Normalized horizontal PV gradient 

EqL - VEC 

deg 

Distance (EqL) from vortex edge center 

EqL - VEI 

deg 

Distance (EqL) from inner vortex edge 

EqL - VEO 

deg 

Distance (EqL) from outer vortex edge 

Single-Level Fields 

Latitude 

deg 

MLS Level 2 reported latitude 

Longitude 

deg 

MLS Level 2 reported longitude 

Dyn Tropopause 

hPa 

Dynamical tropopause pressure 

TG Tropopause 

hPa 

WMO tropopause pressure 


cific method of definition. During spring and fall the auto- 
mated determination also can become problematic. Figure 2 
shows the vortex edge as a function of time, versus sPV and 
ACE CH4, for December through March 2005, in the upper 
and lower stratosphere; in mid-March 2005, the vortex be- 
gan to break up in a major final warming [e.g., Manney et at., 
2006]. In the lower stratosphere, the vortex edge is marked 
by a very well defined region in sPV and CH4 (and other 
trace gases and dynamical markers such as windspeed, PV 
gradient and EqL, not shown) in January through late March, 
becoming slightly less distinct only after mid-March during 
the final warming. Behavior in the middle stratosphere (not 
shown) is similar, with more indication of mixing in CH4 as 
the vortex became more variable in late February. In the up- 
per stratosphere, the region of the vortex edge is much less 
defined throughout the winter because of its complex struc- 
ture. Thus, while the automated vortex edge criterion (and 
any similar criteria) is robust for studies of the winter mid- 
dle and lower stratosphere during relatively quiescent con- 
ditions, for studies involving the upper stratosphere and/or 
fall and spring when the vortex is more variable, due caution 
should be used in applying such automated definitions, and 
their appropriateness should be checked against the physi- 
cal conditions. Note further that the conditions under which 
this, or any other, automated vortex edge definition is robust 
are exactly those under which any reasonable definition of 
the vortex edge (including, e.g., a specific sPV contour) will 
provide a comparable value. 


The WMO (temperature gradient) tropopause is defined 
as the lowest altitude where the temperature lapse rate drops 
below 2 K/km and remains below that for at least 2 km. 
The WMO tropopause is calculated using the algorithm of 
Reichler et al [2003]. For the SO instruments, it is calcu- 
lated from the meteorological analyses’ temperatures after 
they have been interpolated to the SO instrument’s measure- 
ment location and vertical grid; for MLS it is calculated 
before interpolation on the native grid of the meteorolog- 
ical analysis and interpolated to the MLS times/positions. 
The “dynamical” tropopause is defined by the 3.5 x 10“ 6 PV 
contour in the extratropics, joined to the 380 K isentropic 
surface in the tropics or subtropics where that PV contour 
rises above this level. The 3.5 x 10" 6 PV definition is con- 
sistent with the value found to be appropriate for the ex- 
tratropics by Highwood and Berrisford [2000] and Schoe- 
berl [2004]. Tropopause altitude is saved for the SO instru- 
ments, and pressure for MLS, according to the native ver- 
tical grid of the data. Figure 3 shows the WMO and dy- 
namical tropopause altitude for MLS and the five SO in- 
struments during January 2005. The WMO and dynami- 
cal tropopause calculations agree quite well in the summer 
hemisphere (SH) through midlatitudes in the winter hemi- 
sphere (NH); as noted previously [Highwood and Berrisford, 
2000, and references therein], there is often a deep, nearly 
isothermal layer in the polar winter, and thus the WMO 
tropopause is often not very well defined. At the latitudes 
covered by the SO instruments, the range of tropopause val- 
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ues sampled by them is typically similar to that sampled by 
MLS; fewer values from the SO instruments in the low end 
of the range may simply reflect inaccuracy in the conversion 
of the MLS tropopause location from pressure to altitude. 
The SO data also show fewer high values at the WMO trop- 
ical tropopause, which may be related to the finer vertical 
grid used to locate the tropopause for the SO instruments. 

3,2. Sampling Issues 

The DMPs help us to explore the effects of sampling and 
coverage on a variety of comparisons and analyses of the 
satellite data by providing versions of the same fields as sam- 
pled by different instruments. 

Even with the relatively dense coverage of MLS, the sam- 
pling can have a substantial effect on how we view the atmo- 
sphere. Figure 4 shows PV from the 1x1.25° GEOS-4 anal- 
yses in the lower and upper stratosphere, and maps of the 
same fields gridded from the MLS GEOS-4 DMPs. (MLS 
fields are mapped on a 2x5° grid using a weighted average of 
all the points in a day within a specified distance of the grid- 
points, as is commonly done for “quick-look” maps.) The 
day shown is a case where the MLS sampling captures some 
small scale features quite well, e.g., the intrusion into the 
vortex near 30°E at 490 K, and the very narrow double fila- 
ment drawn off the vortex near 120-1 80°E, 30°N at 1700 K. 
However, much of the small-scale structure inside the vortex 
at 1700 K is either distorted or not apparent in the MLS DMP 
fields, and other small-scale features are smeared out (e.g., 
the small high/low PV dipole near 310°E, 30°N at 490 K). 
When comparing features in PV with MLS trace gas obser- 
vations, it is often helpful in assessing agreement to view 
the PV fields as sampled by MLS. With the sparse coverage 
of the SO instruments, sampling effects are always an im- 
portant consideration. Figure 5 shows the latitude sampled 
by each of the five SO instruments used here as a function of 
EqL in the lower stratosphere during the 2004-2005 NH win- 
ter period that will be the focus of time-series comparisons 
shown here, along with the sPV as sampled by those instru- 
ments. As seen in Figure 2, and apparent in the location of 
strong gradients, the vortex edge is near 1.3-1.5xl0~ 4 s -1 
at this level. HALOE and SAGE II sampled very little vor- 
tex air during the winter. SAGE II coverage moves into the 
vortex in late February, largely because of the increasing 
variability in vortex shape and position [e.g., Manney et ai 
2006], but note that, although its measurements in the latter 
half of March reached the highest EqLs, they were taken at 
much lower geographical latitudes than those of the other in- 
struments measuring inside the vortex at this time (^55-65° 
as opposed to ^68-85°). Throughout the winter, the three 
instruments that measured primarily at high latitudes (ACE, 
POAM III, SAGE III) also covered high EqLs by measuring 


at quite different latitudes at a given time. While this should 
not be a large factor in comparing them for species that fol- 
low the vortex closely (e.g., N2O), it can be significant for 
species that may be affected by processes that do not align as 
well with the vortex (e.g., O3, temperature), as will be seen 
in section 4. 

Figure 5 also shows a wide variation in density and com- 
pleteness of coverage of the vortex. This may have a strong 
effect on our perception of vortex-averaged values when 
looking at results from different instruments. The extreme 
example is when comparing MLS (which has nearly com- 
plete coverage of the vortex with dense sampling every day) 
with the SO instruments. Figure 6 shows the average lati- 
tude of all measurements inside the vortex (defined by the 
1.4 xlO” 4 s” 1 sPV contour, c.f., Figure 2) during the NH 
2004-2005 fall and winter, from MLS and ACE measure- 
ments. As expected, MLS average latitudes are typically 
from ~65° to ~72°N during the winter (except at lowest 
levels where vortex is smaller), decreasing dramatically dur- 
ing the vortex breakup starting in early March as the vor- 
tex shrinks and its remnants move away from the pole. In 
contrast, the latitude from the SO instruments is determined 
almost entirely by their measurement patterns, as seen for 
ACE. In section 4, we show examples of vortex-averaged 
trace gas comparisons between MLS and SO instruments, 
using the full MLS sampling and a sampling coincident with 
that from an SO instrument to demonstrate how this may af- 
fect intercomparison results. 

4 Comparisons of MLS with Solar 
Occupation Data Using DMPs 

In the following, we show comparisons of several MLS 
species with SO instrument data using the DMPs to assist 
in comparing measurements in the same air masses. The 
main types of displays are EqL/0 vertical sections, EqL/time 
series, and vortex averages. The EqL/0 plots are produced 
by taking weighted averages around each gridpoint in EqL, 
0 and uncertainty (as provided in each instrument’s data 
files) as described by Manney et al. [1999, 2001]; a Gaus- 
sian weighting function is used. The EqL grid spacing is 
5°, with all points out to four half-widths included to maxi- 
mize coverage for sparse SO instruments; thus, many MLS 
points are included in the average for each gridpoint. The 
0 grid is equally spaced in ln(0), with ~0.09 spacing, cor- 
responding to an ~3-km grid, comparable to the resolution 
of many of the instruments, and generally coarser than the 
data grids; two half- widths are included in 0. EqL/time plots 
are produced in the same way, but with gridding in time in- 
stead of 0. Time gridpoints are 12 UT each day. The SO 
instruments have a 5° latitude half-width and 3-day time 
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half-width, with 3 half-widths included for each. For the 
qualitative time-series comparisons shown here, we have re- 
tained 2.5° latitude, 1.5 day time half- widths for MLS, with 
three half-widths included in each direction. The philoso- 
phy of the uncertainty weighting is similar to that of Man- 
ney et al. [1999]: A relatively narrow half- width is chosen, 
with many half- widths included, so that only points with ex- 
tremely large uncertainties are excluded but the most precise 
data are given much higher weight. 

For vortex averages, the vortex edge is defined by the 
1.4 x 1CT 4 s" 1 sPV contour; very similar results are obtained 
using the EqL from the vortex edge center. For both SO and 
MLS, the measurements are interpolated linearly in ln(0) to 
the plot grid, and all points in the vortex on a day are binned 
and averaged. For averages of MLS measurements coinci- 
dent with SO observations, the criteria are ±1° latitude, ±8° 
longitude, and 12 hours. 

Timeseries plots (EqL/time and vortex averages) are done 
using MLS vl.5 data, to show consistency of time evolu- 
tion and specific features; comparisons that involve only 
ACE and MLS use GEOS-4 DMPs, others use MetO DMPs. 
EqL/0 plots are done using 6- to 8-day periods for which all 
MLS data have been reprocessed in v2.2 for more quantita- 
tive comparisons; these periods include 25-3 1 January 2005, 
and 10-16 March 2005. EqL/0 plots use GEOS-5 DMPs 
for comparisons with ACE, and MetO for all others; where 
the satellite data are recommended for use at those levels, 
MLS/ACE comparisons are shown to the higher altitudes al- 
lowed by the GEOS-5 DMPs. 

4.1. Ozone 

Stratospheric O3 monitoring is one of the primary focuses 
of the Aura MLS experiment. O3 measurements are avail- 
able from MLS and all of the SO instruments considered 
here. Froidevaux et al [2007b, this issue] and Jiang et al 
[2007, this issue] discuss quality and validation of v2.2 MLS 
O3, as well as differences between vl.5 and v2.2. 

Figure 7 shows a comparison of the time evolution of 
O3, using MLS vl.5 data, in the lower stratosphere during 
the 2004-2005 NH winter. Differences between vl.5 and 
v2.2 O3 at this level are small: V2.2 is, on average, ~2-4% 
lower in the lower stratosphere [Froidevaux et al , 2007b, 
this issue]. 2004-2005 was arguably the winter with the 
largest lower stratospheric chemical O3 loss of any NH win- 
ter [WMO, 2007, and references therein]. The morphology 
of O3 was also complex, being highest near the vortex edge 
even before chemical loss began; this complicates separation 
of dynamical and chemical effects and makes assessment 
of O3 loss more problematic [e.g., Manney et al , 2006]. 
It is clear that the O3 evolution in the vortex and along its 
edge agrees quite well, both quantitatively and qualitatively. 


between the four instruments, MLS, ACE, POAM III, and 
SAGE III, with good coverage of the vortex. SAGE II val- 
ues along the vortex edge and in late March also agree quite 
well. POAM III and SAGE III show slightly higher values 
in 40-60°EqL than MLS, SAGE II and HALOE; this may 
be a sampling effect, as the POAM III and SAGE III mea- 
surements are all taken at higher latitudes where high O 3 
values are not completely aligned with the vortex [see, e.g., 
maps in Manney et al, 2006] . ACE vortex interior values are 
slightly higher overall than those from the other instruments, 
but with a very similar amount of O3 decrease between Jan- 
uary and March to that of the other instruments. Exami- 
nation of late-January/mid-March differences as a function 
of EqL and 0 (Not shown) from SO instruments and MLS 
v2.2 data indicate very similar amounts of decrease in the 
lower stratospheric vortex (a signature of chemical loss) for 
each of the instruments, ~0.9 ppmv in MLS, POAM III and 
SAGE III and ^1.1 ppmv in ACE, with maximum changes 
very closely colocated near 70°EqL and 490 K. Detailed cal- 
culations [Singleton et al, 2007] for this winter showed very 
good agreement in O3 loss estimated from these instruments, 
especially when MLS vl.5 data were sampled similarly to 
the SO instruments. 

Figure 8 shows vortex averages of MLS and ACE O 3 , 
along with a vortex average of the MLS points that are coin- 
cident with ACE. This clearly demonstrates some sampling 
effects. In the lower stratosphere, below ~600 K, MLS and 
ACE vortex averages agree very well; these are levels where, 
in absence of heterogeneous chemistry, the lifetime of O 3 is 
very long, and when chemical loss does occur, the patterns 
are generally aligned with the vortex. In the middle and up- 
per stratosphere, in contrast, the full MLS average shows 
considerably lower values than ACE, and significant differ- 
ences in time evolution, especially in the episodes of higher 
vortex O3 (associated with periods of increased wave activ- 
ity). By averaging only those points coincident with ACE 
measurements, most of the qualitative differences in time 
evolution are eliminated, and the MLS values at the O 3 max- 
imum are increased, lessening the difference from ACE. The 
coincident MLS vortex average above ~800 K is typically 
~0.3-0.6 ppmv less than that for ACE, generally consistent 
with the results of Froidevaux et al. [2006] for hemispheric 
averages of MLS vl.5 and ACE v2.1 data. 

Figures 9 and 10 compare v2.2 MLS O 3 with each of 
the SO instruments during March 2005. First, Figure 9 
shows MLS O3 mapped in EqL/0 using the GEOS-5 and 
the MetO DMPs, and the difference between them; this 
demonstrates that the differences in gridding from using dif- 
ferent DMPs are small compared to differences seen below 
between instruments. ACE comparisons in Figure 10 use 
GEOS-5 DMPs and extend to higher levels; other compar- 
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isons use MetO DMPs for MLS as well as SO instruments. 
Figure 10 shows very good agreement between all instru- 
ments in the lower stratosphere and around the tropopause; 
differences are typically less than 0.2 ppmv, and, even at 
these low O3 values, less than ~5%. The exceptions are 
high biases with respect to HALOE and SAGE II in the NH 
tropics and mid-EqLs at the lowest levels, and low biases 
with respect to ACE at the same levels near ±40°EqL, at 
the edges of its coverage. In the middle to upper strato- 
sphere (^ 900-1600 K) near 40~60°EqL, MLS is higher 
than ACE, POAM III and SAGE III (high-latitude-sampling 
instruments) and lower than SAGE II and HALOE (low- 
latitude-sampling instruments). This is the region where 
“low-ozone pockets” (LOPs) [Manney et al. , 1995b; Har- 
vey et al , 2004 ] are common: High O3 from low latitudes is 
drawn into and trapped in the anticyclone for many days. 
When thus confined at high latitudes, O3 relaxes chemi- 
cally toward the lower equilibrium values for these lati- 
tudes (whereas air outside the anticyclones is continually 
being mixed with higher-03 air from low latitudes). The 
pattern seen in MLS data, with highest values apparently 
excluded from this EqL band, is characteristic of the mor- 
phology seen when a LOP is present [e.g., Manney et al , 
1999 , 2001 ]. LOPs were observed continuously from early 
December 2004 through late March 2005 [V. L. Harvey et 
al., “Low-ozone pocket definition using EOS-MLS obser- 
vations”, in preparation, 2007, hereinafter Harvey et al, in 
prep.); these LOPs were sampled by MLS and all SO in- 
struments examined here. Figure 1 1 shows a time-series at 
850 K for MLS vl.5, POAM III and HALOE; MLS vl.5 
and v2.2 typically differ by <2% in the middle stratosphere, 
[Froidevaux et al . , 2007b, this issue]. HALOE provided ex- 
cellent observations of LOPs in January and February (green 
region outside vortex). POAM III sampling of LOPs was 
best in March when the anticyclone was at very high lat- 
itudes. The “pulses” seen in multiple extra-vortex ozone 
minima (especially in MLS where sampling issues are mini- 
mal) are associated with anticyclone merger events [Harvey, 
et al, in prep.). None of the SO instruments sample the 
LOPs and anticyclone region as fully as MLS. Conversely, 
there are cases where an SO instrument samples a tongue of 
very high O3 drawn up between the vortex and anticyclone 
that is missed by MLS [it Harvey, et al., in prep.]. EqL/0 
comparisons like those in Figure 10 constructed using only 
MLS observations coincident with the SO instruments (not 
shown) reduce the differences in the middle stratosphere to 
less than ~ 0.4 ppmv, and in some cases change the sign 
(e.g., MLS is slightly lower than ACE in this region), con- 
clusively demonstrating that these differences result from the 
interplay of sampling, transport and chemistry in this region. 
These differences are of a similar character to those found 


by Manney et al. [2001], and vary with the time period ex- 
amined, depending on the presence and sampling of LOPs 
and low-latitude filaments; similar patterns are seen in EqL/0 
comparisons during the 2005 SH late winter (not shown). 

In the upper stratosphere (and lower mesosphere for 
ACE), Figure 10 shows MLS O3 to be low with respect 
to ACE and SAGE III; this difference is global and persis- 
tent across the time periods examined. Similar differences 
were seen by Froidevaux et al. [2006] with respect to ACE; 
Walker et al. [2005] also noted that ACE vl.O O3 was high 
with respect to other instruments from about 40 to at least 
60 km, and this difference persists in v2.2 ACE data. Froide- 
vaux et al. [2007b, this issue] show a few percent increase in 
v2.2 over vl.5 in the upper stratosphere, with a continuing 
low bias (by ~5-10%) with respect to v2.2 ACE data in v2.2 
MLS data. The high bias of MLS with respect to ACE in the 
lower mesosphere above ^>3000 K is also persistent across 
the periods examined. 

Figure 12 shows O3 at the tropopause from MLS v2.2 
data and each of the SO instruments for 18-day periods in 
January-February and September 2005. Since O3 gradients 
across the tropopause are extremely strong [in fact, O3 is 
sometimes used to define the tropopause, e.g., Bethan et al ., 
1996], this is a very sensitive comparison; also, in this re- 
gion of very low O3 and increasing possibility of clouds, the 
retrievals are difficult. Overall, the values agree remarkably 
well, especially considering that measurements at the same 
latitude from different instruments may have been taken sev- 
eral days apart. In January, all instruments show cases of 
negative O3 values, indicating difficulties in retrieving such 
low O3 values, but most of the values are ~0.1 to 0,4 ppmv 
at the dynamical tropopause, consistent with previous es- 
timates [e.g. Bethan et al., 1996; Pan et al., 1997]. The 
MLS average at the dynamical tropopause, near 0.15 ppmv, 
agrees quite well with the center of the distributions for 
ACE and SAGE HI, and is slightly higher than HALOE and 
POAM III. Ozone varies more at the WMO tropopause, 
and shows a very large scatter in the winter polar regions, 
where that definition of the tropopause is most ambiguous 
(Figure 3); MLS, ACE, and SAGE III all show large (near 
2.0 ppmv) values here in January/February, when the tem- 
peratures and tropopause height were highly variable [e.g. 
Sc hoe her l et al., 2006] and the WMO definition can be am- 
biguous; a similar plot for March 2005 (not shown), when 
the WMO tropopause is better defined, shows much lower, 
more typical, values in the Arctic. A very similar pattern is 
seen in the Antarctic winter, but with less scatter in the O3 
values for all instruments, likely due to more symmetric con- 
ditions in the SH winter resulting in more uniform sampling. 

The preceding discussion of O3 indicates overall very 
good agreement between MLS and all of the SO instruments 
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studied, but highlights many of the sampling issues that are 
important in interpretation of sparse measurements. Map- 
ping in EqL is very valuable in scientific studies with SO 
data, allowing examination of the full range of conditions 
sampled and averaging of air from similar air masses. But, 
particularly with chemically active species such as O 3 , great 
care must be taken to understand how the sampling may af- 
fect the interpretation of the results. Comparisons such as 
those shown above not only help assess the consistency of 
MLS and SO data, but also can assist in understanding and 
interpreting sampling effects to get the maximum informa- 
tion from the sparse SO datasets for scientific studies and for 
validation comparisons with other sparse measurements. 

4.2. Long-lived Tracers 

Three species measured by MLS that are commonly used 
as tracers of transport are N 2 O (in the lower and middle 
stratosphere), H 2 O (throughout the stratosphere and into 
the mesosphere) and CO (mid-upper stratosphere and meso- 
sphere; also upper troposphere and lowermost stratosphere). 
Figure 13 shows vortex averages in the 2004-2005 NH fall 
and winter of these three species, from ACE and MLS. 
The vertical range for CO and H 2 O extends just across the 
stratopause; N 2 O is shown only to the upper stratosphere. In 
general, the strong descent from the mesosphere through the 
middle stratosphere is represented very well in both datasets 
in CO and FLO, with the ACE sampling capturing the extent 
and timing of it. Variations in the lower stratosphere (high- 
lighted in N 2 O) are not quite as fully represented in ACE, 
but after about 15 January, the magnitude of decreases and 
then increases in N 2 O is correctly estimated by ACE. During 
early October through early January (across the gap in ACE 
data), mixing into the vortex resulting from strong wave ac- 
tivity led to increases (decreases) in vortex-averaged N 2 O 
(H 2 0) in the lower through middle stratosphere. When MLS 
is sampled like ACE, however, a significant overestimate of 
these increases results, since ACE sampled near the edge of 
the vortex at this time, where N 2 O (FLO) is higher (lower) 
and mixing in of outside air is greater. “Vortex-averages” 
such as this from SO instruments are often used to estimate 
descent rates for polar processing studies [e.g., Nassar et al . , 
2005] - situations such as this could lead to a substantial 
overestimate of descent from ACE measurements after the 
period of strong mixing. 

Comparing ACE with coincident MLS measurements 
shows very good quantitative agreement in CO and ELO 
in most regions, but slightly lower MLS H 2 O near the 
stratopause. N 2 O shows excellent qualitative agreement in 
time evolution, but a slight overall low bias of MLS with 
respect to ACE, as was noted for NH means of coincident 
profiles by Froidevaux et al. [2006]. V2.2 MLS N 2 O is on 


average about 10% higher than vl.5 between 46 and 2 hPa 
[Lambert et al., 2007, this issue], which would reduce the 
bias seen here. 

4.2.1. N 2 O Figure 14 shows times series as a function 
of EqL at 490 and 850 K of MLS vl .5 and ACE N 2 Q. Very 
good qualitative agreement is seen at both levels, and quan- 
titative agreement at 850 K. At 490 K, MLS is considerably 
lower inside the polar vortex and slightly higher outside, 
with stronger gradients across the vortex edge. Sampling ef- 
fects could be a factor in the differences outside the vortex, 
since ACE is sampling primarily at high latitudes, where the 
ambient N 2 0 is lower, and may miss sampling lower latitude 
air that is not precisely correlated with the PV fields used in 
the EqL mapping. The low MLS bias in the vortex will be 
reduced somewhat in v2.2 [Lambert et al , 2007, this issue], 
and sampling may contribute to it as well (but see below) be- 
cause the sparse sampling of ACE may miss extreme values 
that are represented in the MLS N 2 0 fields; however, a low 
bias in the vortex with respect to ACE appears to be a persis- 
tent feature in both versions of MLS data, and is consistent 
with the results of Lambert et al [2007, this issue]. 

EqL/0 sections comparing MLS v2.2 and ACE data are 
shown in Figure 15 for the NH in March; similar patterns of 
difference appear in plots for the SH winter period examined 
(not shown). These show the persistent low bias of MLS 
in the polar vortex, and in the summer polar lower strato- 
sphere. Similar plots using only the MLS data coincident 
with ACE (not shown) have only slightly smaller biases in 
the lower stratosphere, indicating that sampling differences 
are not the primary factor producing this pattern. Lambert 
et al [2007, this issue] also show a low bias in this region in 
hemispheric averages of differences between MLS and ACE 
coincident profiles. Elsewhere, localized differences are up 
to -15 ppbv, but most areas are within ~5 ppbv, with ACE 
slightly lower than MLS. 

4.2.2. ELO The time evolution of vl.5 MLS H 2 0 is 
compared to ACE in Figure 16 in the lower stratosphere. 
Vortex H 2 0 values show very good agreement, but, once 
again, gradients across the vortex edge are weaker in ACE, 
leading to generally higher midlatitude ACE H 2 0. The con- 
sistency of this with the patterns seen in N 2 0 (Figure 14) 
suggests that the ACE sampling is the primary factor in this 
behavior. In contrast, ACE and MLS H 2 0 agree very well 
inside the vortex, suggesting that the low bias of MLS N 2 0 
with respect to ACE discussed above and seen in Figure 14 
reflects a real difference between the fields being compared. 
In the middle stratosphere (not shown), MLS and ACE H 2 0 
agree well quantitatively throughout the NH winter; in the 
upper stratosphere (not shown), time evolution is similar in 
in ACE and MLS, but ACE has slightly higher overall val- 
ues. POAM III H 2 0 (not shown) is noisier and displays 
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some artifacts; nevertheless, similar time evolution is appar- 
ent, but with values biased high with respect to MLS and 

ACE at all levels. 

Figure 17 shows EqL/0 comparisons of MLS v2.2 H 2 O 
with ACE and HALOE during March 2005. MLS and ACE 
agree extremely well throughout the range, with differences 
typically under 0.2 ppmv (less than 5%) and maximum lo- 
calized differences up to 15%. MLS is higher than HALOE 
by ~0.2-L0 ppmv (~5-25%) throughout the range shown. 
SAGE II and POAM III H 2 0 fields for this period (not 
shown) are somewhat noisy, but similar comparisons with 
them show SAGE II to be ~5-15% lower than MLS be- 
tween 500 and 1600 K, and POAM III to be ~5-15% higher 
than MLS between 400 and 1600 K. Each of these biases is 
consistent with hemispheric biases with respect to vl.5 H 2 0 
data in Froidevaux et al. [2006]; Lambert et al. [2007, this 
issue] show little overall bias between vl.5 and v2.2 H 2 0 in 
the stratosphere. 

4,2.3* CO The chemical lifetime of CO varies from a 
few weeks to many months, making it a useful tracer of 
transport [e.g., Solomon et al., 1985; Allen et al., 1999] in 
many situations; Filipiak et al. [2005] showed initial results 
for MLS CO, and Clerbaux et al. [2005] for ACE CO. As 
shown in Figure 13, starting in fall very high CO descends 
from the mesosphere into the winter stratospheric vortex. 
Figure 18 shows the time evolution of MLS vl.5 and ACE 
CO in the upper stratosphere. Quantitative agreement be- 
tween ACE and MLS is good, as is also the case in the 
middle stratosphere (not shown). The combination of sparse 
sampling and longer time averaging used in the EqL/time 
gridding for ACE are the primary factors responsible for the 
small qualitative differences in short-lived features near the 
vortex edge in late January and early February. Froidevaux 
et al [2006] showed little overall bias between ACE and 
MLS vl.5 CO in the stratosphere, but vl.5 CO often has 
substantial oscillations in the vertical throughout the strato- 
sphere (reduced in v2.2, Pumphrey et al. [2007, this issue]) 
that may result in biases at individual levels. 

Figure 19 shows an EqL/0 comparison of MLS v2.2 and 
ACE CO for March 2005; similar patterns are seen in Jan- 
uary and September (not shown). The MLS CO measure- 
ment is noisy, as seen in the variations in mid-EqLs in the 
upper stratosphere, and v 2.2 still shows some remaining ver- 
tical oscillations in the stratosphere [Pumphrey et al, 2007, 
this issue]; because of this, and the extremely large range of 
CO values, difference plots (not shown) are difficult to in- 
terpret. However, qualitative agreement is good throughout 
the range. Quantitative agreement is best in the Arctic lower 
mesosphere, above about 2000 K (50 km). A high bias of 
MLS with respect to ACE in most regions, but low bias from 
~600 to 800 K, is consistent with profile comparisons shown 


by Pumphrey et al. [2007, this issue]. 

4.3. Other Fields 

4.3.1. HNO 3 HNO 3 is critical in both the activation 
and the deactivation of chlorine, and thus indirectly regu- 
lates ozone destruction in the stratosphere [Solomon, 1999; 
Santee et al, 2004, and references therein]. In the lower 
stratosphere (not shown), time evolution of MLS vl.5 and 
ACE HNO 3 at a function of EqL agrees quite well, but ACE 
misses some of the lowest values in the vortex indicative of 
sequestration in PSCs, as expected for such localized fea- 
tures that are not well-correlated with the vortex. In the mid- 
dle stratosphere (Figure 20, [MLS data are scaled by 0.7 to 
correct an ^30% high bias because of a spectroscopy error in 
vl.5, Santee et al, 2007b, this issue]), qualitative evolution 
appears generally consistent. The low region that develops 
inside the vortex edge in late January through mid-February 
provides an illustrative example of differences related to the 
sparse ACE sampling. Maps of MLS 850 K HNO 3 during 
this period (Figure 21) show the evolution of a partial band 
of low HNO 3 just inside the vortex edge (a similar feature 
appears in v2.2 HNO 3 ). ACE observation locations are over- 
laid, and show that ACE typically measured near the edge of 
this low HNO 3 region where values were slightly higher, and 
the times (e.g., 5 February) when it did sample near the low- 
est values were when that region was farther into the vortex 
interior (thus, at higher EqL); this pattern, and the three-day 
time half- width used for the ACE gridding, result in the faint 
(and distorted toward higher EqL and later time) echo of the 
feature in the ACE data. 

EqL/0 comparisons of MLS v2.2 and ACE HNO 3 in Jan- 
uary 2005 (Figure 22), during a period of strong PSC activity 
[Jin et al , 2006], show good qualitative agreement, includ- 
ing the position and morphology of the primary HNO 3 peak 
and a secondary peak in the polar winter upper stratosphere 
(^900-1300 K) that is a common feature of the HNO 3 distri- 
bution [e.g., Orsolini et al, 2005]. Quantitative agreement 
is fairly good in the levels surrounding the peak (~550 to 
900 K), with MLS typically up to 1.0 ppbv (~5-20%) lower, 
consistent with the results of Santee et al [2007b, this issue]. 
Somewhat larger differences are seen in the mid to low EqL 
lower stratosphere, below ~500 K, and above ~900 K, with 
MLS biased low. Similar patterns are apparent in the other 
periods examined. ACE does not capture the region of low 
HNO 3 in the NH lower stratospheric vortex indicative of se- 
questration in PSCs; this is, again, a sampling effect, since 
the region of PSCs is not aligned with PV contours and is 
limited enough that it is often missed by the ACE sampling. 

4.3.2. HC1 HC1 is a primary reservoir species of Cl 
[e.g., Solomon , 1999]. Froidevaux et al [2007a, this issue] 
discuss detailed validation of MLS v2.2 HC1 and changes 
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from vl.5. EqL/time series evolution of MLS vl.5 and ACE 
HC1 in the lower stratosphere during the very cold 2004- 
2005 NH winter (not shown), when chlorine was activated 
for a prolonged period [e.g., Dufour et al., 2005], indicates 
very good quantitative agreement in the evolution of the re- 
gion of low HC1 associated with chlorine activation, and in 
other features of the time evolution. Figure 23 shows EqL /8 
comparisons of MLS v2.2, ACE, and HALOE HC1 in Jan- 
uary 2005. Agreement between ACE and MLS is very good, 
typically better than 5%, everywhere above ~550 K, consis- 
tent with Froidevaux et al. [2006, 2007a, this issue]. In the 
region of low HC1 in the NH vortex, ACE is higher than 
MLS by up to ~0.4 ppmv, resulting from ACE sampling not 
covering the location of minimum values sampled by MLS. 
HALOE HC1 is ~5- 25% lower than MLS through most of 
the range, consistent with Froidevaux et al. [2006, 2007a, 
this issue]. 

4.3.3. Temperature The use of vortex centered (i.e., 
EqL) coordinates is not an obvious choice for temperature 
from a scientific point of view, since temperature is not, in 
general, expected to be closely spatially correlated with the 
vortex [e.g., Manney et al ., 1996a, 2003; Mann et al., 2002]. 
However, we are often interested in temperatures within the 
vortex, and the comparison between MLS and SO instru- 
ments can help us understand how the sparse sampling of 
the SO instruments may affect studies of processes depend- 
ing on temperature. Also, the EqL coordinate still provides a 
method of comparing SO measurements with other datasets 
over a broader range of conditions than can easily be done 
with individual coincidences or zonal means. The DMPs 
also provide a convenient way of comparing the satellite 
datasets’ temperatures to those in the meteorological anal- 
yses [e.g., Pawson etal, 2007, this issue]. 

Figure 24 shows vortex-averaged MLS vl.5 and ACE 
temperatures, and the average of MLS temperature measure- 
ments coincident with ACE locations. In the middle and 
upper stratosphere, all of these agree quite closely, indicat- 
ing that not only are ACE and MLS temperatures in good 
agreement, but also that the ACE sampling is representa- 
tive of a vortex average. In the lower stratosphere, the full 
MLS vortex average shows more continuous periods of low 
temperatures (as ACE sampling moves in and out of cold 
regions), but higher minimum values (when ACE samples 
cold regions, the “average” is lower than the average includ- 
ing both cold and warmer regions sampled by MLS). ACE 
minimum values are slightly lower than the coincident ones 
for MLS; MLS vl.5 temperatures have a known high bias 
with respect to meteorological analyses and many correla- 
tive datasets [e.g. Froidevaux et al., 2006; Schwartz et al, 
2007, this issue] throughout the stratosphere; Schwartz etal. 
[2007, this issue] show that v2.2 MLS data typically have a 


smaller, low bias with respect to the same datasets. 

Figure 25 compares EqL/0-mapped MLS and ACE tem- 
peratures in September 2005, during the SH late winter. 
ACE is higher than MLS v 2.2 almost everywhere, by a much 
larger amount than the small bias noted by Schwartz et al. 
[2007, this issue]. During this period, there was strong 
wave activity throughout the stratosphere, with the vortex 
distorted and shifted off the pole, and a dipole pattern of 
low/high temperatures across the vortex edge on the side 
with the anticyclone (Figure 26). ACE measurements were 
in the outer part of the vortex, so tended to miss the cold- 
est region and sample near the warmest region on most days 
during the period, resulting in the apparent large high bias 
with respect to MLS. This shows how strong a factor sam- 
pling is in comparisons of temperature, and demonstrates 
the care that must be taken in interpreting temperatures in 
vortex-centered coordinates. 

5. Summary and Outlook 

Derived meteorological products (DMPs) have been cal- 
culated and made available for Aura MLS and the solar oc- 
cultation (SO) instruments ACE, HALOE, POAM II and III, 
and SAGE II and III. These products are fields from and cal- 
culated from gridded meteorological analyses that are inter- 
polated to the observation locations of the satellite instru- 
ments. The DMP fields include 0, horizontal winds, PV, 
horizontal PV gradients, EqL, geopotential height, vortex 
edge criteria, and tropopause locations; DMP files for some 
instruments include additional potentially useful quantities 
such as temperature gradients and line-of-sight PV and tem- 
perature gradients. The DMPs have been calculated from 
the Met Office stratosphere-troposphere assimilation dataset 
for all instruments; in addition, DMPs from GEOS-4 and 
GEOS-5 are available for ACE and MLS (GEOS-4 for MLS 
vl.5 and GEOS-5 for MLS v2.2), and from NCEP/CPC data 
for SAGE II. DMPs are not limited to the Aura mission 
time-period, but have been calculated for the entire mission 
for each instrument. DMPs are provided (see Appendix) 
with the hope that other researchers will find them useful 
both for validation and in science studies combining multi- 
ple datasets. 

Examples comparing MLS and SO data were used to il- 
lustrate sampling issues. Situations where these effects can 
be important include: 

• Comparison of fine-scale transport features (e.g., fila- 
ments and intrusions into the vortex) as represented in 
high-resolution gridded meteorological analyses and 
in MLS data. 

• Studies of O 3 in the middle stratosphere in low ozone 
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pockets, where the combination of chemical and dy- 
namical processes results in O 3 distributions that are 
not closely aligned with the vortex. 

• O 3 distributions in the lower stratosphere when the 
morphology of O 3 is complex and chemical loss is 
taking place. 

• Studies using N 2 O to estimate descent, where the 
movement of SO sampling from vortex edge to vor- 
tex center over time can result in an overestimate of 
descent using “vortex averages”. 

• Localized features, such as PSCs (reflected in depres- 
sions in gas-phase HNO 3 ), that are captured by MLS 
but not by sparse SO sampling. 

• Comparisons during polar winter when temperatures 
are often not concentric with the vortex, such that 
sparse SO sampling misses extreme values. 

SO data are extremely valuable for understanding the at- 
mosphere and monitoring its long-term changes, because of 
their high vertical resolution and precision, the availability 
of long, calibrated datasets, and the potential for detecting 
many species (e.g., ACE). Even greater value is realized in 
studies combining SO and other data, such as UARS MLS 
[e.g. Manney et al ., 1999, 2000; Randel et al. , 1999] and 
Aura MLS [Braathen et al. , 2006; Singleton et al , 2007, 
M. L Santee, et al., “A study of stratospheric chlorine par- 
titioning based on new satellite measurements and model- 
ing”, submitted to Journal of Geophysical Research - At- 
mospheres]. The spatial sparsity of the SO data makes un- 
derstanding sampling effects critical; using the DMPs, the 
sampling limitations have been studied extensively here, to 
help guide interpretation of the SO data and combination of 
them with other datasets. 

Comparisons of MLS and SO data using the DMPs show 
good agreement in time evolution of MLS and SO observa- 
tions in examples from the 2004-2005 Arctic winter; specif- 
ically; 

• Time evolution and values of O 3 in the lower strato- 
spheric vortex agree well between MLS vl.5 and 
ACE, HALOE, POAM III, SAGE II and SAGE III, 
and all instruments show a similar amount of decrease 
indicative of chemical O 3 loss. 

• ACE and MLS vl.5 CO and KLO provide a simi- 
lar picture of the descent of mesospheric air into the 
stratospheric vortex, both in magnitude and timing. 

• When sampled similarly, the time evolution of ACE 
and MLS vl.5 CO, LLO, N 2 O, O 3 and temperature 
averaged in the vortex agree very well. 


• MLS vl.5 N 2 O shows excellent qualitative agreement 
with ACE in EqL/time evolution in the middle and 
lower stratosphere; a low bias of MLS with respect to 
ACE in the lower stratosphere is a persistent feature. 

• MLS vl.5 and ACE H 2 O show excellent agreement 
in EqL/time evolution in the middle and lower strato- 
sphere; likewise for CO in the middle and upper 
stratosphere. 

• MLS vl.5 HNO 3 (scaled to correct a known bias) 
agrees well with ACE HNO 3 in EqL/time evolution, 
with differences explained by sampling effects. 

Comparisons of EqL/0 sections of MLS v2.2 and SO data 
show overall good agreement, including: 

• O 3 agrees to within 5% for for most comparisons, ex- 
cept in regions where differences in MLS and SO sam- 
pling of low O 3 pockets result in larger discrepancies; 
ACE data in the upper stratosphere and lower meso- 
sphere are slightly low with respect to MLS and other 
instruments. 

• MLS v2.2 N 2 O shows a persistent low bias with re- 
spect to ACE in the polar lower stratosphere. Else- 
where, differences are typically within ~5 ppbv, with 
localized differences up to ^15 ppbv. 

• MLS v2.2 and ACE EbO agree very well throughout 
the stratosphere, typically within 5%. HALOE H 2 O 
shows a persistent low bias with respect to MLS of 
r^5-10%, consistent with previous studies. 

• MLS v2.2 and ACE CO compare very well qualita- 
tively, and overall quantitatively; however, MLS CO 
is noisy and still contains some vertical oscillations in 
the stratosphere that compromise detailed quantitative 
agreement. 

• MLS v2.2i and ACE HNO 3 agree very well qualita- 
tively, and to within 5-20% in the region surrounding 
the peak, ~550 to 900 K. At higher and lower levels, 
MLS shows a low bias with respect to ACE. 

® MLS v2.2 and ACE HC1 agree well throughout the 
stratosphere; HALOE HC1 is ~5-25% lower than 
MLS and ACE. 

In addition, comparisons of SO with v2.2 MLS O 3 at the 
tropopause agree well between all of the instruments. These 
results are consistent with more conventional validation re- 
sults presented in other papers in this issue, and help to ex- 
tend those results to a broader range of conditions. 
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Further work to add to and improve DMPs is planned, 
including calculating DMPs for all SO datasets from GEOS- 
4 (and GEOS-5 when that reanalysis is complete) data, 
and calculating DMPs for SAGE I and possibly for other 
datasets, such as UARS MLS and UARS Cryogenic Limb 
Array Etalon Spectrometer. We hope to continue operational 
production of the DMPs for MLS and ACE for the dura- 
tion of those missions. A procedure is available for users to 
request DMPs calculated at locations that they define (see 
Appendix); the system for this is being used to produce 
DMPs for ACE validation campaigns [e.g., Walker et al , 
2005; Kerzenmacher et al , , 2005] conducted in Eureka in 
2004 through 2007 and planned for 2008. We are using the 
DMPs in several studies combining SO datasets and MLS 
data, and hope that they will be useful to other researchers 
in similar efforts. 

Appendix: DMP Access, File Format, and 
Usage 

Table 3 shows the locations for SO DMP access via 
anonymous ftp from mls.jpl.nasa.gov; these can also be ac- 
cessed from the MLS webpage (http://mls.jpl.nasa.gov) un- 
der http://mls.jpl.nasa.gov/research/meteorology.php. The 
ACE DMPs are provided to the ACE Science Team and 
distributed by them. For information on these prod- 
ucts, please contact ACE Mission Scientist Peter Bemath 
(info @ acebox.uwaterloo.ca). 

DMPs can also be obtained by request 

at user-defined (LD) times/locations; see 
http ://mls . j pi .nas a. go v/docs/genericDMPde scription . tx t 
for information on how to make such a request and the 
format for the input files. 

Aura MLS DMPs are publicly available at 

http://mls.jpl.nasa.gov/dmp. Users must first request 
access and agree not to redistribute the MLS DMP files. 
This access request is not intended as a restriction, but rather 
as a means to gauge community uses for these products and 
to provide users with updates and information on changes 
or problems. MLS DMPs can be downloaded for vl.5 and 
v2.2, GEOS-4 or 5 and MetO, for the full Aura mission. 

The formats of the DMP files are designed to follow as 
closely as possible the formats of the datasets they are pro- 
vided for. Information on the file formats and links to sample 
read software are provided on the ftp or web sites along with 
the DMPs. 
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HALOE 
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Relative to pub/outgoing/manney after logging in via FTP 
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3 <y>=4~digit year, <d>=3-digit day-of-year 
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Figure 1. Monthly averages of equivalent latitude 
(EqL)/potential temperature (0)-mapped (left to right) wind- 
speed (ms -1 ), normalized (see text) horizontal PV gradient, 
and distance in EqL from the vortex edge center (°EqL). 
EqL/0 mapping is from derived meteorological products 
(DMPs) calculated from GEOS-4 data for MLS locations, 
during (top) January 2005 and (bottom) March 2005. White 
space in EqL from vortex edge plots indicates that the vortex 
is not defined anytime during the month in these regions. 




Equivalent Latitude / deg 


Figure 2. Timeseries of distance in EqL from vortex edge 
center (colors, °EqL) as a function of (left) sPV (10 -4 s -1 ) 
and (right) ACE CH 4 (ppmv) at (top) 1700 K and (bottom) 
520 K, for ACE measurement locations during December 
through March 2005. Vortex edge location values are from 
GEOS-4 DMPs. 


This preprint was prepared with AGU’s macros v5.01, with the 
extension package ‘AGU ++ ’ by P. W. Daly, version 1.6b from 1999/08/19. 
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Figure 3. (Left) WMO and (right) dynamical (see text) 
tropopause altitudes from MetO DMPs for MLS (range 
shown as grey shading, white line average), ACE (green), 
HALOE (red), SAGE II (blue), SAGE III (purple) and 
POAM III (gold) for all measurements in January 2005. 
Tropopause pressures from MLS DMPs are converted to al- 
titude using a scale height of 7.0 km. 
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Figure 5. EqL/time plots of latitude (°, left) and sPV 
(10~ 4 s -1 ) sampled by (top to bottom) ACE, POAM III, 
SAGE III, HALOE, and SAGE II during December 2004 
through January 2005 EqL on the 490 K isentropic surface. 


Figure 4. PV (1(T 4 K m 2 kg" 1 s" 1 ) maps from GEOS-4 
analyses (top row) and gridded (see text) from the MLS 
GEOS-4 DMP PV field (bottom row), at 490 K (left) and 
1700 K (right). Projection is orthographic, with 0° longitude 
at the bottom; latitude range is from 0° to 90°N. 
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Figure 6. Vortex-averaged (within the 1.4 xlO -4 s -1 sPV 
contour from GEOS-4 DMPs) latitude (°) sampled by (top) 
MLS and (bottom) ACE from 15 October 2004 through 
31 March 2005. 



MANNEY ET AL. : Meteorological Products for Aura Intercomparisons 


22 




40 1 


30 i 

o 

g. 

20< 

50 

40 

30 

20 


1400 


800 


Figure 8. Vortex-averaged (within 1.4xl0 -4 s _1 sPV con- 
tour from GEOS-4 DMPs) O 3 (ppmv) for 3 October 2004 
through 31 March 2005 from (top) MLS vl.5, (center) MLS 
vl.5 measurements coincident with ACE (see text), and (bot- 
tom) ACE. Vertical range is 400 to 2500 K. 
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Figure 10. Global EqL/0 cross-sections for 10- 
16 March 2005 comparing MLS v2.2 O 3 with (top to bot- 
tom) ACE, POAM HI, SAGE III, HALOE, and SAGE II. 
Left panels show the SO instruments’ O 3 and right panels the 
MLS-SO instrument O 3 difference (ppmv). For MLS/ACE 
comparisons, GEOS-5 DMPs are used for EqL/0 mapping; 
others comparisons use MetO DMPs. MLS/ACE compar- 
isons are shown from 360 through 4430 K; others are 360 
through 2500 K (the top limit of the the MetO DMPs). 
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Figure 7. EqL/time series for 1 December 2004 through 3 1 March 2005 of O 3 (ppmv) at 490 K from (left to right, top to 
bottom) MLS (vl.5), ACE, POAM III, SAGE III, HALOE, and SAGE II. EqL mapping is done with MetO DMPs. Overlaid 
contours are sPV (10 -4 s -1 ), with three contours in the vortex edge region shown as solid lines. 
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Figure 9. Global EqL/0 cross-sections for 10-16 March 2005 comparing MLS v2.2 O 3 (ppmv) mapped using GEOS-5 (left) 
and MetO (center) DMPs, and the difference (ppmv, right). Vertical range is 360 through 4300 K. 
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Figure 11. EqL/time series for 1 December 2004 through 
3 1 March 2005 of O 3 (ppmv) at 850 K from (top to bottom) 
MLS (vl.5), POAM III, and HALOE. Mapping is done us- 
ing MetO DMPs. Overlaid contours are sPV (10 -4 s _1 ), 
with three contours in the vortex edge region shown as solid 
lines. 
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Figure 12. O 3 (ppmv) at the (left) WMO and (right) dynam- 
ical tropopause for MLS v2.2 (range shown as grey shading, 
average white line), ACE (green), HALOE (red), SAGE II 
(blue), SAGE El (purple) and POAM III (gold) for all mea- 
surements during (top) 24 January-10 February 2005 and 
(bottom) 10-27 September 2005. Tropopause locations for 
interpolation of O 3 are from MetO DMPs. 
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Figure 13. As in Figure 8, but for (top to bottom) CO (ppbv), H 2 0 (ppmv), and N 2 0 (ppbv). Vertical range for CO and HiO 
is 400-2500 K; for N 2 0, 400-1600 K. Vortex edge sPV is from GEOS-4 DMPs. 
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Figure 16. As in Figure 7, but for MLS vl.5 (top) and ACE 
(bottom) EbO (ppmv) at 490 K, using GEOS-4 DMPs. 
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Figure 14. As in Figure 7, but for MLS vl.5 (top) and ACE (bottom) N 2 O (ppbv) at (left) 490 and (right) 850 K, using 
GEOS-4 DMPs. 



Figure 15. Global EqL/0 cross-sections for 10-16 March 2005 comparing MLS v2.2 with ACE N 2 O Left panel shows MLS 
v2.2 N 2 O, center panel, ACE N 2 O, and right panel the MLS-ACE N 2 O difference (ppbv). GEOS-5 DMPs are used for EqL/0 
mapping. Vertical range is 400 through 1800 K. 
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Figure 18. As in Figure 7, but for MLS vl.5 (top) and ACE 
(bottom) CO (ppbv) at 1700 K, using GEOS-4 DMPs. 
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Figure 19. As in Figure 15 (with difference panel omitted), 
but for MLS v2.2 CO (ppbv) compared with ACE, using 
GEOS-5 DMPs. Vertical range is 36CM-300 K. CO contour 
interval is on a log scale, from 20 to 10,000 ppbv. 
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Figure 17. As in Figure 15, but for MLS v2.2 H 2 0 (ppmv) compared with (top) ACE and (bottom) HALOE. GEOS-5 DMPs 
are used for ACE comparison and MetO for HALOE comparison. Vertical range for ACE comparison is 360-4300 K, for 
HALOE comparison, 360-2500 K. 
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Figure 20. As in Figure 1 1, but for MLS vl.5 (top) and ACE 
(bottom) HNO 3 (ppbv) at 850 K, using GEOS-4 DMPs. 



Figure 21. Maps of MLS vl.5 HNO 3 (ppbv) at 850 K for 
28 and 30 January and 2 and 5 February 2005, with ACE 
observation locations (white dots) overlaid. PV contours in 
the vortex edge region are overlaid in black. Projection is 
orthographic, with 0 ° longitude at the bottom; latitude range 
is from 0 to 90°N. 
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Figure 22. As in Figure 1 5 , but for MLS v2.2 and ACE HNO 3 (ppbv) for 25-3 1 January 2005, using GEOS-5 DMPs. Vertical 
range is 400-1600 K. 
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Figure 23. As in Figure 15, but for MLS v2.2, ACE, and HALOE HCI (ppbv) for 25-31 January 2005. ACE comparison 
uses GEOS-5 DMPs, HALOE comparison uses MetO DMPs. Vertical range is 400-2500 K. 
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Figure 24. As in Figure 8, but for MLS vl.5 temperature 
(K). Vertical range is 400-2500 K. GEOS-4 DMPs are used 
to define the vortex edge. 
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Figure 26. Maps of MLS v2.2 Temperature (colors, K) on 
22 September 2005 at (left) 490 K and (right) 850 K. Over- 
laid contours are 1.4 and 1.8 xlO -4 s -1 . White dots show 
ACE observation locations. 
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Figure 25. As in Figure 15, but for 17-24 September 2005 MLS v2.2 and ACE temperatures, using GEOS-5 DMPs. Vertical 
range is 360 to 4300 K. 


